Criterion on remote clocks synchronization within a Heisenberg scaling accuracy 
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We propose a quantum method to judge whether two spatially separated clocks have been syn- 
chronized within a specific accuracy a. If the measurement result of the experiment is obviously a 
nonzero value, the time difference between two clocks is smaller than a; otherwise the difference is 
beyond a. On sharing the 2iV-qubit bipartite maximally entangled state in this scheme, the accuracy 
of judgement can be enhanced to a ~ n/(uj(N + 1)). This criterion is consistent with Heisenberg 
scaling that can be considered as beating standard quantum limit, moreover, the unbiased estimation 
condition is not necessary. 
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Introduction.- Quantum entanglement, a distinctive 
feature of quantum mechanics, is at the heart of appli- 
cations in distributed systems, e.g., quantum key distri- 
bution and clock synchronization. Clock synchronization 
with high precision is a fundamental and an important 
problem in that it is crucial for many modern technolo- 
gies and researches, such as global positioning system 
(GPS), long baseline interferometry, synchronous data 
transfer, gravitational wave observation (LIGO), tests of 
theory of general relativity, and distributed computation. 
There are two standard methods for synchronizing two 
spatially separated clocks in the frame of special the- 
ory of relativity. One is based on Einstein's synchro- 
nization procedure which uses an operational linc-of-sight 
exchange of light pulses between two spatially separated 
clocks The other one is based on the internal time 
evolution of quantum systems, like Eddington's infinitcs- 
imally slow clock transport |2|. The quantum clock syn- 
chronization method based on the strength of sharing 
prior entanglement has also been proposed in 0] , and has 
been generalized to several multiparty clock synchroniza- 
tion protocols [4-6]. Independent of the parties's knowl- 
edge of their relative locations or of the properties of the 
intervening media, these procedures utilize the instanta- 
neity of wavefunction collapsing after the measurement is 
performed on the shared entangled states. Since the pro- 
cess of distributing entanglement is adiabatic, these pro- 
tocols arc tantamount to Eddington protocol. There are 
some relevant experimental works done on the quantum 
clock synchronization protocols: an experiment focusing 
on the quantum clock synchronization implementation 
has also been reported in Ref. 0; in addition the pro- 
gressive techniques of multi-photon entanglement Q are 
basic and promising for the realization of these protocols. 

Quantum entanglement-enhanced parameter estima- 
tion that plays a vital role in quantum metrology uses 
quantum mechanical property to enhance the sensitiv- 
ity of the measurement of classical quantities. It has 
been pointed that the standard quantum limit 1/y/N, 
where TV is the number of particles used in the measure- 



ment, can be beaten by using the coherent light with 
squeezed vacuum In the study of quantum metrol- 
ogy, quantum Fisher information theory and quantum 
Cramer-Rao bound based on the statistical distance of 
states have been proposed and developed in EMI- The 
NOON state has been demonstrated to be able to achieve 
a phase sensitivity saturating the Heisenberg limit 1/N 
[16| . Some related strategies have been proposed to per- 
form a high precision in quantum metrology framework 
whilst man y expe riments have also been per- 
formed on this topic 2(J- HI ■ One direct and natural idea 
to apply this technique is in the clock synchronization. 
Chuang [2(| has presented a high efficiency quantum tick- 
ing qubits handshake protocol which allows two remote 
clocks to be synchronized independent of message trans- 
port time; and a similar protocol has been proposed to 
beat the standard quantum limit (27j . 

In this Letter, we relate the quantum clock synchro- 
nization protocol to the problem of estimating an un- 
known parameter. We investigate the performance of 
the bipartite maximally entangled spin-zero singlet in 
the scheme of two-clock synchronization and offer a stan- 
dard to judge whether two spatially separated clocks have 
been synchronized in a specific accuracy. This criterion 
is practical and is consistent to Heisenberg scaling, addi- 
tionally it does not rely on the unbiased estimation con- 
dition which is a fundamental hypotheses in the quantum 
Fisher information theory 

General framework of quantum clock synchronization. - 
Suppose two spatially separated parties, Alice and Bob, 
rested on the same reference frame, both possess high- 
precision clocks, such as Cs atomic clocks, running at 
exactly the same rate. They do not agree on a common 
time at the same readout, for example twelve o'clock. 
The difference of time Y between their clocks can be 
expressed as 



Y = t B -t A \ 



Alice and Bob have the same readouts- 



(1) 



In a quantum scheme, in order to eliminate the rela- 
tive phase that may emerge during qubits transport- 
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ing to the spatially separated locations, the entangled 
states should be distributed to Alice and Bob adiabat- 
ically. After the entanglement distribution, Alice and 
Bob respectively perform measurements on all of their 
qubits simultaneously when their clocks point to the 
same readout. We choose a normalized entangled state 
|V') = X)i jPij\iA)\3B) which is merely changed with an 
overall unobservable phase under the unitary evolution 



results is 



u AB {t) 



-iH A t/H 



-iH B t/H 



\%a) and \jb) are or- 



thonormal basis of measurements which satisfy the com- 
pleteness Ylij Na)|.7,b)0b|(*a| = I where I denotes to 
the identity. Suppose Alice performs the measurement 
before Bob and obtains a result \ia) with probability 
P(ia) = \Pij\ 2 , then the collapsed state evolves as 



(2) 



Then Bob will obtain 



where U kj = (k B \e- tH ^ h \j B 

the result \k B ) with probability P{k B \iA 
Thus, Fisher information theory and Cremer-Rao bound 
can be utilized in this clock synchronization situa- 
tion while the estimation is asymptotically unbiased. By 
comparing the ratio of observed measurement outcomes 
with probability distribution that is determined by the 
parameter Y, two issues may prevent one estimating Y 
with a high precision. First, the number of experimental 
trials is finite, so the ratio of measured outcomes may 
deviate from the distribution; only when the number v 
is large enough, could the estimation be unbiased and 
the Creamer-Rao bound be reached. Second, a one-to- 
one mapping P(^\Y) -h- Y between the probability dis- 
tributions and parameter is essential. In this Letter, we 
assume that Alice and Bob respectively perform measure- 
ments X = 1 0) (0 1 — 1 1) ( 1 1 on all of their own qubits simul- 
taneously when their own clock points to a specific value, 
where \x) = ^ E^oC" 1 )""*!*) ( x = 0,1), and |0),|1) 
are the orthogonal eigenstates of each qubit, which have 
the identical Hamiltonian H satisfying H\Q) = Eo\Q), 
H\l) = Ei\l), and u = (Ex - E )/h > 0. One way to 
implement these ticking qubits in experiment is to put 
some spin- 1/2 particles into the magnetic fields with the 
same field strength. 

Quantum clock synchronization with Bell state and 
GHZ state.- At first, we suppose that Alice and Bob 
share Nv pairs entangled qubits with form |^( - )) 
: 1 1 10) ) / a/2 Q which is invariant under unitary evo- 
lution e - tH At/h ^ e -iH B t/h _ Alice and Bob perform mea- 
surements expressed as an operator f = X(tA) ® X{t B ) 
in the Heisenberg picture, and it can be described as a set 
of positive operator valued measurements with element 
E(Q = I^K^I ® \y)(y\ while £ = (x, y) where x,y = 0,1. 
Then, the probability distribution of all the measurement 



P(£\Y) - Tr[E(£)p(Y)] = - (4,^ cos , 



' ^sin 2 



and the Fisher information is calculated as 
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Ty =Y J P W) 



dlnP(£\Y) 



dY 
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(3) 



(4) 



where p(Y) is the density matrix of pure state e 



,-t.ffr 



e .„ d (t+r)/ft|^(-)^ p = u\Y\, and 5 x , y is Kronecker's 
delta. In addition the average of the measurement opera- 
tor can be calculated as (/) y = g(£)P(j£\Y) = — cos /3 
where = (— l) x+y . If \Y\ < tt/lj holds, we can obtain 
the difference of time \Y est \ from the observed expecta- 
tion value / = j^jj a ft er measurements and 
classical communication. The sign of Y est can be deter- 
mined by the outcomes of Alice's and Bob's measure- 
ment because the one who firstly performed Nv times 
measurement would get a duel results with probabili- 
ties P(|0)) = P(|T)) = 1/2. Furthermore, the uncer- 
tainty of estimation of Y could reach the Cramer-Rao 
bound SY est = 1/ (cjy/ NvFy) = 1/(uj\/Nv) which is the 
standard quantum limit. Therefore, in this scheme Al- 
ice and Bob can synchronize their clocks with accuracy 
l/{uy/Nv). 

Some researches have also focused on using the quan- 
tum entanglement strategies and employing GHZ like 
states, |GHZ) = (|0)^|l)| w + |l>? w |0>f w Uv^, to en- 
hance the precision of parameter estimation [l6Hl9( • It is 
easy to verify that GHZ state is merely changed with an 
overall phase under the unitary evolution (e~' lHt / h )® N ® 
-im/h^®N _ rpj^ pj-Qbabii^y distribution in this proto- 



(e 

col takes the form: P(£'|F) = ^r[l + g'(C) cos(A/3)], 
where g'(C) = { — l)^^ Xk+Vk ' 1 while the symbol is 
£ = (xi,- ■- ,x N ,yx,-- ■ ,y N ), x l ,y j = 0,1; and /' = 
X(tA)® N <8> X(t B )® N is the measurement operator in the 
Heisenberg picture. The average of the operator is cal- 
culated as cos(N/3), and Fisher information Ty is N 2 uj 2 . 
Considering that the probability distributions and expec- 
tation value are all function with periodicity 2n/N, thus 
one can unambiguously obtain |y es t| from the observed 
expectation value after measurements and classical com- 
munication only when the condition \ Y\ < tv/(uiN) is sat- 
isfied. The sign of Y est can also be determined by the out- 
comes of Alice's and Bob's measurement because the one 
who firstly performed Nv times measurement would get 
the probability PflxT- ■ ■ x£)) = 1/2 N and P(\xi)) = 1/2. 
When v is large enough, the uncertainty can attain 
the Cramer-Rao bound SY est = \j\JvTy = \/{ujNy/u) 
which has a Heisenberg scaling accuracy. Despite of this 
optimal local distinguishability in the Hilbert space, GHZ 
states are inappropriate to obtain more advantageous in- 
formation from any values of parameter Y in this single 
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procedure since the condition \Y\ < tt/(wN) is required 



Quantum clock synchronization with bipartite maxi- 
mally entangled states. — We next consider a scheme 
which exploit different entanglement resource. The bi- 
partite maximally entangled spin-zero singlet has been 
proposed as a resource for quantum-enhanced metrology 
}29j , with the following form: 



\x) 



V2J 



L= y (-iy- M 



\J,M) Z! a\J,-M) ZjB (5) 



M--J 



where J = N/2, and \J,M) Z is a completely symmet- 
ric normalized state (Dicke state) with (J — M) qubits 
being |0) and (J + M) qubits being There is an 

explicit mapping between the two-symmetric entangled 
state and the direct product of N maximally entangled 
states, which is presented in [30[, then one obtains 



2 N/2 



VN + 1 

2 N/2 



Nly/N- 



AiB ai ■ ■ 



AnB 



,(6) 



where I is the identity operator on Hilbert space TL = 
{|0),|1)}, S = T,m=-j\J>M) z (J,M\ is the symmetric 



projector that maps states in onto its symmetric 

subspace Hf 1 *, a denotes to a permutation. After the 
adiabatic distribution of Nu pairs entanglement l^ - )), 
Alice or Bob can perform the symmetric projector S to 
obtain v pairs \x). Because |#(->) = (|01) - \10))/y/2 
changes only with an overall unobscrvablc phase under 
any unitary evolution of form U <8> U in two-qubit space, 
then this singlet has the rotational invariancc property 
under unitary evolution U® N <E> U® N and has identical 
expression in any spin basis, e.g. zh>ih>j; when ignor- 
ing the overall phase. Another proof has been presented 
in Ref . (3ll | , and this invariance property has been tested 
in the experiment j32[. These bipartite maximally en- 
tangled states play an important role in the quantum 
information distribution and concentration [33J435|. Re- 
cently these states used in our scheme have been gen- 
erated experimentally by using stimulated parametric 
down-conversion and have been used in the 1 to 3 + 2 in- 
formation distribution [H, |3(| . Additionally, some other 
experiments also produce such entanglement and real- 
ize the quantum information distribution [371 ] . Next, we 
show that these technologies can also be utilized to im- 
plement our scheme for quantum clock synchronization. 
The pure state evolved as 

^ e -iH A t A /h^®N ^ r e -iS B tB/h\ig>N \j\ 
= (H® N <g> M® N )(I® N <g> ( e -<VH fl HflHB//l)®Arj| x j 

= (M® N ® M® N )(I® N ® U B (n/2, (3, -n/2f N )\ X )(7) 




FIG. 1: The functional relation between the expectation value 
f(/3) and the parameter /?. 



1 1 

1 -1 



is 



where the overall phase is ignored; Ma.b = ^ 

the Hadmard matrix and the unitary operator U(a, /3, 7) 
is expressed by three Euler angles in the basis {|0), |1)} 
as the following: 



£/(«,/?, 7) 



:|p»(«+7)/2 



£p-i(a-7)/2 



. sm | e ^-7)/2 cos | e - t (a+ 7 )/2 



(8) 



Moreover, according to group theory of the irreducible 
representation, we obtain an analytical expression of the 
unitary operator in 7V-qubit space [HI 



U (a, /?, 7 r v | J M) 



M' 



< Ma + M '^d J M , M (l3)\JM').(9) 



Thus, we can obtain the probability distribution of mea- 
surement outcomes £' = (xi,--- ,xn ,yi, ■ ■ ■ ,tjn), with 
Xi,yj = 0, 1; 



P(£'|Y) 



{2J+\)C J 2 J M C J 2 J M 



(10) 



where (J — M) is the number of in {xi,--- ,xn} 
while (J — M') is the number of in {yi, ■ ■ ■ , yjv}- As 
before the expectation value of measurement operator 
/' = X(t A )® N ® X(t B )® N is calculated as 



fifi) ■■= (f'h 



E (-D 

M,M'=-J 

(-l) N sin(A + l)/3 



N+M+M' 



2J+ 1 



N + 1 



sin ft 



(11) 



The function f((3) against its argument (3 is shown in 
Fig. [1] with different numbers of qubits. Furthermore, 
Fisher information reads 



T Y = 



j 

E 



dY M',—M 



(13)] AJ(J+1)uj 2 



RLM' = -J 



2J+1 



(12) 
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with which it is straightforward to find that lower bound 
5Y est = ^j3/{ujyjN(N + 2)v) obviously breaks the quan- 
tum standard limit and performs a Hcisenbcrg scaling 
accuracy. 

Nevertheless, this scheme has some special proper- 
ties that differ from the previous schemes. Consider- 
ing that /(/?) is clearly "peaked" around (3 — with 
width ~ n/(N + 1) (sec Fig. [J), we can confirm that the 
uncertainty of Y est could reach the Cramer-Rao bound 
1/ ' [uiy/v^Fy) for a large number v when the condition 
\Y\ < n/((N + 1)oj) is satisfied. Although the bipar- 
tite maximally entangled spin-zero singlet fails to gain 
more advantageous information of the parameter Y from 
the observed expectation value against the GHZ state, 
more importantly we will acquire a quantum enhanced 
criterion to judge whether two remote clocks have been 
synchronized with the accuracy tt/ (oj(N +1)) even when 
the number v is not large enough, i.e., the unbiased es- 
timation hypothesis is not fulfilled. When the expec- 
tation value / = i Yh=i flte) obtained from the out- 
comes of experiments after classical communication sat- 
isfies |7- {-l) N \ < 1- 1/V2, one obtains 
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Prob \Y\ < 



(N + l)u) 



> Prob 



f-(fh 



< 



> 1 - 2e" 



1 

(13) 



which is Chernoff bound [39| ■ For example, suppose that 
|7 - (-1)^1 < 1 - 1/V2 and v = 10, we can infer that 
the inequality \Y\ < tt/((N + 1)lu) holds with fiducial 
probability larger than 84%. 

Conclusions. — In many cases of quantum metrol- 
ogy, quantum Cramer-Rao bound can only be achieved 
asymptotically ,i.e., it holds for unbiased estimation with 
infinite number v and zero error SY es t — > |4oT ]: and 
unfortunately this problem also exists in some quantum 
clock synchronization strategies. However, applying the 
bipartite maximally entangled spin-zero singlet one can 
obtain a standard to judge whether two spatially sep- 
arated clocks have been synchronized within a specific 
uncertainty even when the number v is not large enough. 
If so, we can step further to obtain the difference between 
two clocks with a Heisenberg scaling accuracy in accor- 
dance to the expectation of the measurement results by 
increasing the number v. 

In conclusion, we propose a novel quantum scheme for 
remote clocks synchronization within a specific accuracy. 
This bound of accuracy scales as the Heisenberg limit 
which is the ultimate limit of precision measurements 
under all conditions. With developments in creating ex- 
perimentally the entanglement resources, this quantum 
scheme of remote clocks synchronization may be imple- 
mented which may possess unprecedent precision. 

We would like to thank Guo-Yong Xiang, Xiang-Ru 
Xiao and Li Jing for useful discussions. This work 
was supported by 973 program (2010CB922904), NSFC 



[0 
[10 

[11 

[12 



[1.3 

[14 

[15 
[16 

[1- 

[18 

[19 

[20 

[21 
[22 
[23 

[24 
[25 
[26 
[27 

[28 

[29 

[30 

[31 
[32 



Electronic address: muliangzhu@pku.edu.cn 

Electronic address: hfan@iphy.ac.cn 

A. Einstein, Ann. Phys. 17, 891 (1905). 

A. S. Eddington, The Mathematical Theory of Relativity, 

2nd edition, (Cambridge University Press, Cambridge, 

1924). 

R. Jozsa et al, Phys. Rev. Lett. 85, 2010 (2000); E. A. 
Burt et al., Phys. Rev. Lett. 87, 129801(2001); R. Jozsa 
et al, Phys. Rev. Lett. 87, 129802 (2001). 
M. Krco, and P. Paul, Phys. Rev. A 66, 024305 (2002). 
R. Ben-Av, and I. Exman, Phys. Rev. A 84, 014301 

(2011) . 

C. Ren, and H. F. Hofmann, Phys, Rev. A 86, 014301 

(2012) . 

A. Valencia et al., Appl. Phys. Lett. 85, 2655 (2004). 

Z. Zhao et al, Science 430, 54 (2004); X. C. Yao et al., 

Nat. Photon. 6, 225 (2012). 

C. M. Caves, Phys. Rev. D 23, 1693 (1981). 

R. A. Fisher, Proc. Camb. Soc. 22, 700 (1925). 

H. Cramer, Mathematical Methods of Statistics (Prince- 
ton University, Princeton, NJ, 1946). 

A. S. Holevo, Probabitistic and Statistica Aspects of 
Quantum Theory (North-Holland, Amsterdam, 1982), 
especially Chaps. III. 2 and VI. 2. 

S. L. Braunstein, J. Phys. A 25, 3813 (1992); S. L. Braun- 
stein, Phys. Rev. Lett. 69, 3598 (1992). 
S. L. Braunstein, and C. M. Caves, Phys. Rev. Lett. 72, 
3439 (1994). 

B. M. Escher et al, Nat. Phys. 7, 406 (2011). 

J. J. Bollinger, W. M. Itano, D. J. Wineland and D. J. 
Heinzen Phys. Rev. A 54, R4649 (1996). 
V. Giovannetti, S. Lloyd, and L. Maccone , Science 306, 
1330 (2004). 

V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. 
Lett. 96, 010401 (2006). 

V. Giovannetti, S. Lloyd, and L. Maccone, Nat. Photon. 
5, 222 (2011). 

T. Nagata, R. Okamoto, J. L. OBrien, K. Sasaki, and S. 

Takeuchi, Science 316, 726 (2007). 

B. L. Higgins et al, Nature. 450, 393 (2007). 

M. Kacprowicz et al, Nat. Photon. 4, 357 (2010). 

I. Afek, O. Ambar, and Y. Silberberg, Science 328, 879 
(2010). 

G. Y. Xiang et al, Nat. Photon. 5, 43 (2011). 

H. Yonezawa et al, Science 337, 1514 (2012). 

I. L. Chuang, Phys. Rev. Lett. 85, 2006 (2000). 

M. de Burgh and S. D. Bartlett, Phys. Rev. A 72, 042301 
(2005). 

G. A. Durkin, and J. P. Dowling, Phys. Rev. Lett. 99, 
070801 (2007). 

H. Cable, and G. A. Durkin, Phys. Rev. Lett. 105, 
013603 (2010). 

Y. N. Wang, H. D. Shi, Z. X. Xiong, L. Jing, X. J. Ren, 
L. Z. Mu, and H. Fan, Phys. Rev. A 84, 034302 (2011). 
J. Schliemann, Phys. Rev. A 72, 012307 (2005). 
M. Radmark, Marek Zukowski, and M. Bourennane, 



5 



Phys. Rev. Lett. 103, 150501 (2009). 
[33] M. Murao, D. Jonathan, M. B. Plenio, and V. Vedral, 

Phys. Rev. A 59, 156 (1999). 
[34] M. Murao, M. B. Plenio, and V. Vedral, Phys. Rev. A 

61, 032311 (2000). 
[35] Y. L. Zhang et at, Phys. Rev. A 87, 022302 (2013). 
[36] M. Radmark, M. Zukowski and M. Bourennane, New J. 

Phys. 11, 103016 (2009); M. Radmark, M. Wiesniak, 

M. Zukowski, and M. Bourennane , Phys. Rev. A 80, 

040302(R) (2009). 
[37] A. Lamas-Linares, J. C. Howell, and D. Bouwmeester, 

Nature (London) 412, 887 (2001); F. Ciccarello, M. Pa- 

ternostro, S. Bose, D. E. Browne, G. M. Palma, and M. 

Zarcone, Phys. Rev. A 82, 030302(R) (2010); H. Yu, Y. 



Luo, and W. Yao, Phys. Rev. A 84, 032337 (2011); A. 
Chiuri, C. Greganti, M. Paternostro, G. Vallone, and P. 
Mataloni, Phys. Rev. Lett. 109, 173604 (2012). 
[38] Z. Ma, Group Theory for Physicists (World Scientific, 
Singapore, 2007); 

fm_y ( 1 ^ [(J+M')!(J-M')!(J+A/)!(J-A/)!] 1 / 2 

U M',M\I J ) — Z^,A 1 ) (J+M'-v)\{J-M-vy.v\(v+M-M'y. A 

(cos §) 2J + M '- M - 2 »(sm f )2-M'+Af_ 
[39] A. S. Nelsen, and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

Cambridge, England, 2000). 
[40] M. Tsang, Phys. Rev. Lett. 108, 230401 (2012); D. W. 

Berry et al, Phys. Rev. A 86, 053813 (2012). 



